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Abstract 

We analyze the behaviour of a system of particles living on a periodic crystal in the presence of 
a magnetic field B. This can be done by involving a periodic potential U(x) and the Rashba inter- 
action of coupling constant k so . By resorting the corresponding spectrum, we explicitly determine 
the band structures and the Bloch spinors. These allow us to discuss the system symmetries in 
terms of the polarizations where they are shown to be broken. The dynamical spin will be studied 
by calculating different quantities. In the limits: k so and U(x) — 0, we analyze again the system 
by deriving different results. Considering the strong B case, we obtain an interesting result that is 
the conservation of the polarizations. Analyzing the critical point Afc i(T = ^\J\ 5 we show that the 
Hilbert space associated to the spectrum in z-direction has a zero mode energy similar to that of 
massless Dirac fermions in graphene. Finally, we give the resulting energy spectrum when B = 
and U(x) is arbitrary. 
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1 Introduction 



Unlike the quantum Hall effect (QHE) pQ that can be resulted from the impurities [2] or the Coulomb 
intercation [3], the spin Hall effect (SHE) [3] is essentially due to something else. In fact, this phe- 
nomena can be seen as consequence of the Rashba or Dresselhaus interactions [3 [6j [7] . They depend 
what kind of spin-orbit coupling is involved and more precisely they are originated from two different 
effects. For the first one is caused by the asymmetry of the confining quantum well but the second is 
a bulk effect [8]. The spin-orbit coupling leads to a spin accumulation in the direction transverse to 
the applied electric field. These show the importance of such interaction and motivate their analysis 
in different contexts. 

Many theories appeared as good candidates to describe different issues related to the Rashba 
interaction. Different questions have been addressed and lead to interesting developments on the 
subject. In particular an analytical and numerical results for the effect of such spin-orbit coupling on 
band structure, transport, and interaction effects in quantum wires when the spin precession length 
is comparable to the wire width have been reported in [9]. The Harper-Hofstadter problem for two- 
dimensional electron gas with Rashba interaction subject to periodic potential and perpendicular 
magnetic field has been studied analytically and numerically [10J. One may cite also other interesting 
papers related to the subject such as [HI [T2J [131 Ej. 

Very recently, an analysis has been appeared dealing with other issues. This concerns the in- 
vestigations of the basic features of a two-dimensional electron gas with Rashba interaction and two 
in-plane potentials superimposed along directions perpendicular to each other [15]. The first of these 
potentials is assumed to be a general periodic potential while the second one is totally arbitrary. A 
general form for Bloch's amplitude is found and an eigenvalue problem for the band structure of the 
system is derived. Moreover the general result is applied to the two particular cases in which either 
the second potential represents a harmonic in-plane confinement or it is zero. It is found that for a 
harmonic confinement regions of the Brillouin zone with high polarizations are associated with the 
ones of large group velocity. 

Motivated by the investigations of different groups cited above and in particular [15], we develop 
our proposal to deal with different issues. For this, we consider a system of particles living on the 
plane (x, z) in the presence of an external magnetic field. To do our task, we include a period potential 
U{x) as well as the Rashba interaction. By splitting the corresponding Hamiltonian, we diagonalize 
different parts to end up with the total eigenvalues and Bloch spinors for the two first transverse 
modes, i.e n = 0, 1. These will be used to show that there is a hidden symmetry in the spectrum. 

On the other hand, by analyzing the polarizations we conclude that such symmetries are not 
preserved for the present system. Moreover, we study the dynamical spin by evaluating the velocity 
and spin current operator, which are used to give the corresponding Hall conductivity. Moreover, 
we discuss two limits where the first is when the Rashba coupling constant and this leads to end up 
with the Landau problem submitted to the periodicity constraint in x-direction. The second one is 
U{x) = that allows us to get the same Landau problem but amended with the Rashba interaction. 

To deeply investigate the basic features of the present system, we treat an interesting case. Indeed, 
for a strong B we analyze the corresponding spectrum and symmetries to end up with different 
conclusions. As an interesting consequence, we show that the polarization symmetries are preserved 
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in such case and this generalize the obtained result in [TS] to the magnetic field. More importantly, 
such polarizations are obtained to be -B-dependent and therefore it can be controlled by different 
values taken by the magnetic field. 

On the other hand, we consider the critical point Afc j(7 = iy^ that appears as a singularity of 
the spectrum to underline its influence on the obtained results. More precisely, by resorting different 
quantities we notice the spectrum corresponding to z-direction is sharing some commun features with 
the graphene. In particular, a zero mode energy will come up in similar way to massless Dirac fermions. 
In the end, we analyze the periodic structures when B is switched off and by considering two limits 
we recover some known results like the Rashba dispersion relation for instance. 

The present paper is organized as follows. In section 2, for later convenience we study the Bloch the- 
ory for one-dimensional system and establish the Hamiltonian for two-dimensional where the Rashba 
interaction is included. This will be generalized by introducing a magnetic field B perpendicular to 
the periodic crystal in section 3 in order to create a confinement in one-direction. Moreover, we split 
the corresponding Hamiltonian into two parts to get an easy diagonalization. In section 4, we treat 
different parts to obtain the eigenvalues and eigenstates and for the total spectrum we restrict our- 
selves to the two first modes. In section 5, we consider two different limits: k so = and U(x) = 
to characterize the system behaviour in such cases. We analyze the polarization symmetries by eval- 
uating different components in section 6 and show that only that in z-direction is broken where the 
remaining are nulls. We discuss the spin dynamical properties by evaluating some relevant quantities 
in section 7. We consider the case when B is strong to deduce different conclusions in section 8. The 
critical point will be analyzed in section 9. We consider the periodic structure of the system when B 
is switched of in section 9. Finally, we conclude our work. 



2 Bloch theory for one-dimensional system 

We start by showing the spinorial structure influence on one-dimensional periodic crystal. This can 
be done by making use of a Bloch theory and characterize the system behaviour. In doing so, let us 
consider a simple Hamiltonian, which is given by 

where U (x) is a periodic potential, such as 

U (x + L) = U (x) (2) 

and L is its period along the x-direction. Clearly, without U(x) we simply have a plane wave in 
one-dimension as solution of the system. 

Note that, H^ D is a simple problem that can easily be solved. Indeed, the eigenvalues e\°^(k B ) and 
eigenstates \l, k B ,cr) are solutions of the equation 

H™\l,k B ,a) =e^(k B )\l,k B ,a). (3) 

They are labeled by the Bloch's quasi-momentum k B that is a discrete value in the Brillouin zone and 
the band index I. The energies ef^kg) are degenerates with respect to the spin index because one 
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can obtain the relations 

^i(M = 4^i(M = 4 u; (^) ( 4 ) 

where e, (/sb) read as 

^ - ^- (5) 

To get the corresponding Bloch wavefunctions, we simply project the states \l, fcs, <r) on the coordinate 
representation {|x, a}}. Doing this process to get 

(x, a'\l, k B ,a) = -^=e ikBX Ul , kB , a (x, a') (6) 

where ui t k B ,a(x, f') is the Bloch amplitude and Lq is the system size. Moreover, we can easily establish 
the relation 

Ul,k B ,a( X ' a ') = Sv',v u i,k B ( x )- ( 7 ) 
Note that, u^ B (x) is also a periodic function and satisfies 

ui,k B {x) = ui,k B (x + L). (8) 

In conclusion, the spinorial structure of the Bloch theory is trivial in the present case. This suggests to 
see what happens if some interaction in terms of spin are included. This issue has been addressed for 
a system in two-dimensions and submitted to a confining potential [15]. This work will be generalized 
to get others interesting results and offer different discussions. 

Before proceding our generalization by including the magnetic field B, it is relevant to fix the 
starting point that is B = 0. This will allow us to underline the difference between these two cases 
and therefore make some comments. In fact, a Bloch theory in two-dimensional periodic crystal with 
in the spin-orbit coupling described by the Hamiltonian 

Hq = w _ _ + + 

2m m 

was considered in the reference |15j where the second term is the Rashba interaction and the confining 
potential is fixed to be 

V{z) = (10) 

Some interesting results have been derived and in particular the polarization symmetries have been 
discussed. In fact, we will show that these results can be recovered as a particular case of the present 
proposal when B is strong and then it can be used to controle the polarization behaviours. 



3 Gauge field coupling 

To deal with different issues we introduce an external magnetic field B to end up with a confining 
potential in similar way to (|10p but with a frequency -B-dependent. This will allow us to obtain an 
appropriate Hamiltonian that will be used to do our task and in particular to generalize the result 
obtained in [15] to the case of B and offer some comments. 
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Let us start by considering a periodic crystal parametrized by (x, z) in the presence of a perpen- 
dicular B. The Hamiltonian for a single particle can be written as 

vf 2 h 2 K 



H 



where the conjugate momentum is 



2m m 



{oxK z - a z ir x ) + U(x) 



- e 7 
ir = p A. 

c 



(11) 



(12) 



Unlike §§§ the Rashba term is actually ^-dependent and there is no external confinement analogue to 
the potential V(z). However, we will see that an analogue form to V{z) can spontaneously be created 
thanks to the emergence of B. Thus, one can expect that this effect will make difference with respect 
to the results obtained in [15j . The corresponding quasi-momentum becomes 



K = k-^-A 

he 



(13) 



where A is the vector potential. These materials will be the starting point to establish a tool in order 
to tackle different problems in this work. 

In this level, it is necessary to fix some quantities involved in the above Hamiltonian. In doing so, 
we choose the Landau gauge 

A = B{z,0). (14) 



This allows (1111) to take another form, such as 



H 



2m 



eB 

Pt+ [Px~ —Z 



h 2 k v 



m 



(JxPz -a z \p. 



eB 



U(x). 



(15) 



One can use different methods to solve this problem. In our analysis, we propose to split H into two 
parts and each one will be treated separately. They are 



where the Hamiltonian is given by 

H {1) = ^( k * + cr z k so ) 2 + U(x) 
and is fixed to be 



h 2 k 2 h 2 k 2 1 feB\ 



2m 



+ 



+ 



— 



H (2) 



2m 2m \ c J 

H k xn 

~&xkz 



heB 



mc 



rn 



(16) 



(k x + a z k so ) z (17) 



(18) 



Clearly, due to the presence of B, the Bloch spinors will be changed from |Z,/ce,cr) to the states 
\l,kB,rj). In fact, they verify the eigenvalue equation 



H\l,k B ,v) = £l,v{ k B)\l,k B ,v) 



(19) 



where the parameter rj will be fixed later. Then, the next step is to solve such equation in order to 
derive the whole spectrum that can be used to deal with different issues and in particular discuss the 
symmetries. 
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4 Eigenvalue problems 



As it will be clear soon, the above splitting is useful is sense that the whole spectrum can be obtained 
in the easiest way. More precisely, we resort different spectrum entering in the game by an immediate 
derivation from the Bloch theory and making an algebraic approach to overcome some difficulties. 
In doing so, we split H^ 1 ' into two forms as well and analyze each part to get its spectrum. As far 
as the whole spectrum is concerned, we restrict ourselves to two first modes those allow us to make 
comparison with the results obtained in [15] . On the other hand, this restriction is argued by dealing 
with the Rashba interaction effect on the band structures [9]. 

4.1 spectrum 

By separating the Hamiltonian H^' into two parts, we show that the first spectrum can be easily 
obtained from that corresponding to the zero magnetic field case seen before. However, the second 
one can be deduced from one-dimensional harmonic oscillator. This can be done by distinguishing 
two different cases: spin down and up. But, it is enough to handle only one of them and other can be 
obtained in similar way to reach the conclusion. 

The simplest way to derive the corresponding spectrum is to divide itself into parts. The 
advantage of this is to refer to each part independently and also to make diffrent comparisons with 
other cases. Then, the first one reads as 

HP = ^{k x + a z k so ) 2 + U{x) - (20) 
2m 2m 

and the second is given by 

„ m h 2 k 2 z 1 (eB\ 2 , _ , , 

H^ = — - + — — z 2 - fuv c (k x + a z k so )z 21 
2m 2m \ c J 

where u> c = — is the cyclotron frequency. It is now clear that the second term in (|21|) plies the role 
of the confining potential (|10p . This will allow us to make contact with one-dimensional harmonic 
oscillator in order to get the spectrum of (|2ip . 



The H x 1] form is easy to handle where its eigenvalues and eigenstates can be derived in similar 
way to those obtained in section 2. In fact, one can find the spectrum 

= 4 0) (kB + ak so ) - ^ (22) 
where a = ±1 is the spin index. It can also be written as 

h 2 

I* = ^ [k% + 2<rk so ) . (23) 
Clearly in the absence of the Rashba coupling, i.e. k so = 0, we recover that corresponds to 



spinless particle. 
In deriving t: 

that is easy to diagonalize. We proceed by distinguishing two cases: spin up and down, i.e. a = ±1, 



In deriving the spectrum, we use some changes to end up with the harmonic oscillator form 
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where we treat only one of them and the second will follow immediately. Then, we start by defining 
new variables for a = 1, such as 



Pz 



^J2mhw c 



Q 



muj c 
~2h 



■z. 



They verify the commutation relation 



[Q,P] 



In terms of P and Q, \ a= i takes the form 



HP\ a=1 = Hu c 



P 2 + Q 2 



2h 

mLO c 



(k B + k so ) Q 



(24) 
(25) 

(26) 



which is sharing some commun features with the harmonic oscillator. To clarify this point, let us make 
an algebraic analysis by introducing the annihilation and creation operators. They can be realized as 



a = Q + iP + Xk,+i, a) = Q - iP + \k,+i 
where the quantity \k,+i is a function of the Bloch and rashba vectors, such as 



A 



n 



2moj c 



(k B + k so ) 



a, a 



It is easy to check the relation 

With these operators, H^\ a= i can be mapped as 

HP\ a=1 = huj, 



ot a + i(i-2A| )+1 ) 



(27) 

(28) 
(29) 

(30) 



To make contact with the harmonic oscillator spectrum, firstly we define a new frequency in terms of 
lo c and Xk,+i- This is 

w +1 = u c (1 - 2\% +1 ) (31) 



which allows us to write (1301) as 



Hjp\ ff=1 = tkj 



+i 



1 " 2A 2 fe)+1 



+ 2 • 



(32) 



Note that, hereafter we assume that the condition Aj^+i ^ ±y \ is fulfilled. In fact, we return to 
analyze this critical point separately and underline its influences on the system behaviour. Secondly, 
we can consider new operators in terms of a and a T involved before. They are 

1 .* 1 



1 " 2A|, +1 



1 " 2X l + i 



They lead to the Hamiltonian 



(33) 



(34) 
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This is nothing but the harmonic oscillator Hamiltonian and therefore the spectrum can easily be 
obtained. Indeed, the eigenvalues are given by 

4 1 2,+l(^)U = ^+l ( n +l) > n = 0,1,2... (35) 
and the eigenstates read as 

In) = ^ |0>. (36) 



To complete the above analysis we need to consider the second case where a = — 1. Indeed, the 
same calculations remain valid and the only thing that should be changed is \k,+i- Its analogue can 
be defined in similar way, such as 



**.-! = "V 2^** -*"°>- (3?) 
It is easy to observe \k,+i and A^-i verify the symmetry relation 

^k,+i(k B ) = -Afc (38) 

Therefore, from the above analysis we end up with the eigenvalues for spin down case, such as 

^l-l( k B)\z = f^-l(n + ^\, n = 0,1,2-.. (39) 

where u+% and u—x satisfy a similar relation to (|38j) . These last spectrum's will be used to analyze 
the critical point. This will allow to make contact with the Dirac fermions in the presence of magnetic 
field. 

At this level, we have all ingredients to get the spectrum. Indeed, we just combine all together 
to obtain those verify the eigenvalue equation 

H^\l,k B ,n,a) = e^(k B )\l,k B ,n,a) (40) 

which is showing that the eigenvalues of are of the form 

= £< i 0) (** + - + ^ ( n + 2 1 ) ■ ( 41 ) 

It is clear that the corresponding eigenstates can be obtained as tensor product, namely 

\l,k B ,n, a) = \l,k B ,a) (g> \n). (42) 

The eigenfunctions can be obtained by projecting the above states on the coordinate represen- 
tation. Thus, they can be found in terms of the Bloch amplitude as 

(x,n',a'\l,k B ,n,a) = —=d n > jn S a > ,ae ikBX % kB +ak s „(x)- (43) 



With this we complete the derivation of the eigenvalues and the eigenstates problems of These 
will be used to derive other quantities and discuss different issues. 
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4.2 Total spectrum 

After getting the spectrum we need to determine that corresponding to in order to end up 
with the total eigenvalues and eigenstates. In this case, the Bloch amplitude acquires a new spinoral 
structure, such as 

{x, n, a\l, k B ,v) = -7=e ikBX u kkB ^{x, n, a) (44) 

where ui t k B ,rj(x,n, a) is also a periodic function, namely 

ui,k B , v (x, n, cr) = ui ikB>r) (x + L, n, a). (45) 

Before proceeding, let us denote by Oi,k B ,r] (n,a) the Bloch spinors in the {l,k B ,n,a} representation, 
where r\ will be fixed later. They can be written as 

Ql,k B ,v ( n ' a ) = C> n ' ks > ( 46 ) 
which allowing us to derive the relation 

(I', k' B ,n, a\l, k B ,v) = 5 V,i^k' w k B & i,k B ,n («, o) . (47) 

To explicitly determine Qiksw (ti, cr), we have to solve the eigenvalue equation fll9|) where still only 
to be treated. In doing so, let us project this equation on the states \l, fee, n, a) to get the Bloch 
spinors as 

[(l',k' B ,n,a\H^\l,kB,n,a'} + (/', k' B , n, a\H^\l, k B , ri, a')] 6^ (n',c/) = 

n',a' 

= £ i,v(kB)h' ,rfk' B ,k B 9i,k B ,r] (n, a) . (48) 

This is actually involving two kinds of matrix elements where each one can separately be evaluated 
because they are related to and H^ 2 \ Indeed, for the first it is easy to obtain 

(l',k' B ,n,a\HW\l,k B ,n',a'} = e^Jk^^^y. (49) 

The second kind can be written as 

ft?k 

(l',k' B ,n,a\H^\l,k B ,n',a') = —{a\a x \a'){n\k z \n')5 V i8u k . (50) 

m 

This contains also two separated matrix elements. As concerning the first one, its is easy to end up 
with the result 

(a\a x \a') = 1 - 5 a>a >. (51) 

Combining all together to write the final equation for the eigenvalues si^ks) and eigenspinors 

8l,k B , v (n,a) as 

4ll,a( k B)$l',lh B ,k B Sa,a' ~ (l ~ S<r,(r>) (n\k z \n ') 8 V jd 'y g >k f 

n ,<r 

= £ irtkB)di',idk> B ,k B 6i,k B , v (n,a) . (52) 

Actually, this is depending to different quantum numbers and can be solved by making use of some 
restrictions. This will be the subject of the remaining subsection. 
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To solve the above equation, we may focus on some cases those allow us to simplify different 
quantities. Indeed, assuming that we have the same band indices I' = I and Bloch's quasi-momentum 
k B = k' B , then ([52]) can be reduced to 

h?k so 



\ ^n,n'^'£^l a (k B ) 



2m 



(l - 5 a y) (n\k z \n') >0i >kBiV (n' ,a') = ei )r] (k B )9i )kB ,r,(n,cr). (53) 



To go further, we should evaluate the second term in (|53p . This can be done by returning to the 
annihilation and creation operators in order to obtain 

2ih 



A-A^ 

This is leading to write k z as 

k z = i(A-A^ 
Now it is easy to get the matrix element of k z as 

(n\k z \n) = ±z<5 n>n /±i« 
Therefore, the eigenvalue equation takes the form 



1-2A2 



(54) 



2h 



(55) 



' muj a 
~2h 



1 1 

n + - =F - 

2 2 



(56) 



2m 



(k 2 B + 2ok so ) + hw a ( n + 



h 2 k so 

T i — (l — ^<T,cr') <>n,n'±l' 



I muja- I 11 
-2fT { n+ 2 T 2 



xOi >kB ,r,(n', a') = ei tV (k B )0i,k B , v (n, a). 



(57) 



Many discussions can be reported on this spectrum as well as its corresponding eigenfunctions. In 
the next, we treat a special case that allows us to make contact with other results and give different 
comments. More precisely, since we are wondering to make comparison with the proposal developed 
in |15| . we restrict ourselves only on the two first transverse modes, i.e. n = 0, 1. On the other hand, 
such approximation has been motivated by analyzing the spin accumulation in quantum wires with 
strong Rashba coupling [9j. This restriction will allow us to simply <\57h and derive interesting results. 

Focussing on n = 0, 1 and a = ±1, (|57p reduces to a diagonalization of 4 x 4 matrix, which is a 
task that can analytically be done. This process leads to fix r\ as a number running on (1, • • • , 4) and 
to end up with two relations between different energies. They are 



^,77=1,2(^5) = sf{k B ) ~ 61,2 (Xb)> e/,r;=3,4(feB) = Ej 4 ' (k B ) + £,l 2A (k B ) 

where e^(k B ) reads as 



e l [H B ) = 
and the quantity £^ 2 {k B ) is given by 



ef\k B + k so )+ef\k B -k so ) 



2u2 



hLJ c . 1 - 2Xi 



h z k 



2m 



612 ( k 



B 



e7(k B ) T^(l~2Xl a 



(58) 



(59) 



(60) 
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The energy e, (k B ) is 



and £ takes the form 



£ i ( k s) = 2 { e i 0) + fc so) - e| 0) (k B - fc so ) j 



VS ( 1_2A t)- 



(61) 



(62) 



We show that there is an hidden symmetry that encoded in the spectrum. Indeed, from the Bloch 
theory it is easy to verify the relation 



ef\k B )=ef\-k B ) 
which allows us to show that the different energies are connected via 

£l, v =i(k B ) = £i, v =2(-k B ), e^ v=3 (k B ) = e l>v=i (-k B ). 



(63) 



(64) 



The corresponding normalized eigenspinors Oi^ B ,r)( n , c) can be obtained by returning to the eigenvalue 
equation. This leads to get 



Ol,k B ,V=lA ~ N l,k B ,r)=l,S,k B ,V=lAi 0l,k B ,Tl=2,3 ~ N l,k B , v =2,3 e i,k B , v =2,3 

where the functions 6^,fc s ,r;=i,4 take the form 




1 



(65) 



'l,k B ,r)=l,A 



\ 



(66) 



and the remaining quantities read as 

/ 



6l,k B ,r)=2,3 — 



\ 





er(M + ¥( 1 - 2A L)±&(^ 
i 





(67) 



Note that, we have involved such notation 



01 



k B ,V 



( 0, )fcB) „(n = O,(7 = +l) \ 
Ol,k B , v (n = 0,cr = -1) 
®l,k B , v (n = l,a= +1) 

V Ol,k B ,r,( n = 1,0- = -1) J 



(68) 



in deriving the eigenspinors and the normalized constats can be obtained as usual. Thus, it is not 
hard to show that there are 

i +i , 2 



N l,k B ,v = J2 p,k B ,r,(n,a) 

n=0 a=-l 

Finally, using the obtained spectrum it is easy to verify the relation 



(69) 



Nl,k B , V =l,4 - Nz,-fc fl ,7j=2,3 



(70) 
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between different sub-bands entering in the game. 

At this level, let us give a comment about the band structures. Indeed, calculating the derivative 
of £j{k B ) with respect to k B at the zero point, we obtain the group velocity 

de 7( k B) { (0),, \ 

-^-\Ms=o=v\ (k so ) (71) 

which is nothing but that corresponds to one-dimensional system. Clearly, starting from the obtained 
spectrum one can observe that the band splitting near k B = is a linear function in terms of the 
Bloch vector k B . It is similar to what has experimentally been observed for two-dimensional electron 
gazes |21j . Our conclusion is in agreement with what has been reached in |15j . 

To close this section let us note that when the magnetic field is switched off, we recover Bloch 
Hamiltonian amended by a free particle living on the z-direction. To get the analysis of [15] simply we 
can add a harmonic potential in the same direction. Our study shows that this potential can simply 
be generated from the confinement due to the presence of B. On the other hand, we will use the 
established tools to discuss the symmetry in terms the polarizations. 



5 Polarization symmetries 

We start by noting that it has been shown that the system described by the Hamiltonian ([9]) is 
preserving the polarization symmetries [15]. Thus, it is natural to ask about such symmetries for the 
present system in order to underline the similarities and differences. To reply this question, we start 
discussing the polarizations by adopting the definition 

P^(kB) = (l,k B ,vWi\hk B ,v) (72) 

where Oi are the Pauli spin operators, with i = x,y,z and r] = 1, 2, 3, 4. Therefore, we need to evaluate 
different quantities entering in the game. 

We start by analyzing each component separately. In doing so, we write the identity operator in 
the {|Z, k B ,n, a)} basis and take into account the structure of the obtained Bloch spinors. This process 
leads to show that its expression in the x-direction is given by 

i +i 

P tv = E E °tk B , v ^ *0 (1 - <W") 0l,k B)V (^ °") = (73) 

n=0 a',a"=-l 

as well as for the y-direction 

i +i 

P S = E E e ik B ,r,(n, a') f" (1 - 6 a , >a „) e l>kBjr ,(n, a") = (74) 

n=0 <j',ct"=-1 

for all I and k B in the Brillouin zone. Clearly, both of polarizations longitudinal and perpendicular 
are nulls, which are in agreement with the analysis of |15j . 

Let us now analyze the polarization along the z-direction. Indeed, for such case we show that (|72p 
reduces to the quantity 

i +i 

p tt = E E °) ° i*b,v&> °)- ( 75 ) 

n=0 cr=-l 
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To explicitly evaluate this form, we need to separate two cases those governed by the four Bloch 
sub-bands and use of course the obtained spinors. Thus, for r\ = 1,4 we find 



pO) _ jvr-l 



and for r] = 2, 3 we end up with 



pi') - N -i 

r l,r) ~ iv Z,fc B ,)7=2,3 



er(*B) + ^(i-2A2)T&(fcfl) 



i 



n 2 



(76) 



(77) 



This shows this component of polarizations has a finite value, which makes difference with respect to 
those seen above. Moreover, using (|70p and the relations 

= -ef (-A*), 6i, 2 (^) / ^(-A*) (78) 
we show that the symmetry relation for the polarizations verifies 



(79) 



This tells us that the polarization symmetries are not preserved for the present system. These have 
been broken by the external magnetic field that interact with the particles in terms of the gauge 
field. Its preservation can be recovered by switching off B and adding a confining potential along 
z-direction [15]. However, as we will in the next, we can have such preservation of the symmetries by 
considering a limit of strong B. 



6 Dynamical spin analysis 

In the beginning, let us note that the present system is sharing some commun features with that 
describing the spin Hall effect (SHE). In fact, to talk about SHE one has to consider the Rashba or 
Dressellhaus spin-orbit coupling as an important ingredient in the Hamiltonian system. This suggests 
to make a dynamical spin analysis of our system in order to underline its behaviour. This can be 
done by determining the velocity and spin current operators, which can be used to write down the 
corresponding spin Hall conductivity <r^. 

To do our task we start by evaluating the velocity operators in different directions. Indeed, from 
the Hamiltonian formalism one can easily obtain the quantities 

T£ X fo^SO ' Z fo^SO /or\\ 

v x = 1 a z , v z = 1 a x (80) 

mm mm 

which will be used to derive the spin current operator. To do this, let us recall that due to the fact 
that spin is not a conserved quantity in the presence of the spin-orbit coupling, there are different 
definitions of such operator [16J. But we adopt that widely used where the spin current operator is 
defined as the difference of the conserved densities for carriers with opposite spins. This is 

h 

Jz v = -r{v z (Xy + o y v z ) (81) 
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leading to the form 

J? = ^*y. (82) 

2m 

Note that, Jz v can be also obtained by using another method. The advantage of this latter is to end 
up with an appropriate form for u^, which we will do soon. 

Reccal that SHE is a manifestation of a collective motion of the spins of particles in an external 
electric field. This phenomenon is characterized by cr^ that is related to that for the quantum Hall 
effect |17| . From linear response formalism, it can be written in terms of the frequency uj as [16j 



jS y _ eh {a\Jy z \a'){a'\vx\a') f(e a ) - /(eg) 



11 

a.a' 



where we have set \a) = \l,kB,n,a) and e a = £/ )J? (A;b). It can be simplified by introducing Jy z in 
terms of the Heisenberg equation, namely 

J? = - l l [H,a z ). (84) 
The corresponding matrix element in the basis \a) can be evaluated as 

(a\Jz y \a) = - j(^a ~ e a /)(a\a z \a). (85) 

Injecting this in ([83]) . we end up with 

s eh s-^ (a\a z \a')(a'\v x \a') 

4L Z e Q — e a / + akj + ^U 

a, a' 

To explicitly determine such conductivity one may use the self-consistent Born approximation in 
similar way to [16]. On the other hand, (|86p can be analyzed by adopting the high and low frequencies 
regimes to reach different conclusions. 



7 Analysis of two cases 

We are wondering to analyze different limits in order to characterize the system behaviour. From 
this section we start by considering the cases where the Rashba coupling constant and the periodic 
potential are nulls. These will allow us to end up with a Landau problem amended either with a 
periodic potential or the Rashba term in the presence of the magnetic field. 

7.1 Without Rashba coupling 

The first one concerns the Rashba coupling constant when is absent. In fact, we return back to the 
former studies and show its restrictions to such case. In doing so, we derive the corresponding spectrum 
as well as discuss the polarizations symmetries. Before doing our job, let us fix the Hamiltonian that 
describes the present limit. This can easily be obtained from (|15|) to get 



1 

2m 



eB 

Px —Z 



+ U(x) 



(87) 
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which is nothing but the Landau problem written in the Landau gage and submitted to the constraint 
of a periodic potential. Clearly, without B we have a Bloch theory along x-direction plus a wave plane 
along z-direction. 

To get the eigenvalues and the eigenstates, we algebraically diagonalize the above Hamiltonian. 
To proceed, let us define the annihilation and creation operators as 



Ib 

V2h 



Pz + l\ Px 



eB 



Ib 

V2h 



eB 

Pz + t\ Px —Z 



where l B = y jrg is the magnetic length. It is easy to verify the commutation relation 



c, c 1 



I. 



In terms of these operators, ()8Tf) can be written as 



ri\ kso=0 = hu c [Jc+-) +u{x). 



Therefore the corresponding spectrum is given by 

1" 



E n \ kso=0 (k B ) = hw c [n + 



\n,k B ), 



n = 0,1,2, 



where the states explicitly read as 



|0), (x\k B ) 



1 



* ikBX u kB {x). 



(89) 



(90) 



(91) 



(92) 



To close this part, we note that the eigenvalues depend only on n and are completely independent of 
the Bloch momentum hk B . Due to the lack of dependence of the energy on k B , the degeneracy of each 
level is enormous. 



7.2 Without periodic potential 

Let us study another case that corresponds to U(x) equal to zero. This will allows us to characterize 
the nature of the present system when U(x) is switched off. In doing so, we give the corresponding 
Hamiltonian, which can be deduced from (|15p as 

1 h 2 k 
H\u{x)=q = 7T- K + T*) Kvr 2 - a z ir x ) . 93 

Before proceeding we note that an analogue operator has been considered in different contexts but 
when only the x an y coordinates are involved, for instance one may see [18J . Here is different because 
of the last term which is including diagonal and off-diagonal matrices. 

To get the spectrum of (j93|) we adopt an analysis similar to that has been reported in subsection 
(4.1) except that we take U(x) = 0. We also note that for each state we can still choose k x freely, 
thus we have degenerate energy levels. If we assume our space is finite, for definiteness we choose 
a rectangle of dimensions L x x L y , the number of states in a level is finite and can be calculated. 
The eigenvalues of p x are quantized as k x = j^-j where j is a quantum number. Therefore, after 
diagonalizing (|93p we end up with the relations 

^=1,2 (i) = £/ + C0 - 6i, 2 0')> £/,r,=3,4(j) = efij) + 62,1 0') ( 94 ) 
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where s^(j) takes the form 



2m \ L y J 



+ ^ c (l-2A 2 CT ) 



and 2 (j) reads as 

The energy ej~(j) is given by 



*r 0) = 



«r(i) T ^(i-»^) 

ft 2 /2vrj 



m \ L 



and the quantity £ is 



where the parameter A,- a becomes 



7i 2 £w, fw, 



3^ 



h [2irj 



2muj c \ L y 



+ crk s 



These lead to derive the normalized eigenspinors 9ij tV (n,a) as 

_i _ _j_ 

Ol,j,TI=lA = N l,j%=l,A,3,n=lAi dl,k B , V =2,3 = N lj* v =2,3 lJ,V=2,3 

where the spinors 0zj i7? =i,4 read as 





1 



and 0^=2,3 are given by 



e r(j) + ^(i-2A 2 (7 )±6 2 (j)' 



V 



i 





(95) 

(96) 

(97) 
(98) 

(99) 

(100) 



(101) 



(102) 



It is natural ask about the polarization symmetries in the present case. The answer can be obtained 
by evaluating different components and to obtain for 77 = 1,4 the result 



Huj, 



er(i)-^(i-2Ay ±&(i) 



and for 77 = 2, 3 we have 
P W _ 

ly l,j,r,=2,3 

By showing the relation 

>(*) 



(103) 



(104) 



^M(7)b(x)=o / -^S=2 >s (-J)b(x)=0- (105) 
we conclude that such symmetries are also not preserved in the present case. This close the analysis 
of the Rashba and periodic potential terms. In the next, we will study to other cases. 
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8 Strong magnetic field limit 



It is interesting to analyze the case where the champ magnetic is strong. Such interest comes from 
the fact that this limit is important for the quantum Hall effect. Thus, it will be a good task to deal 
with such limit and derive different results. 

Recall that B is included in the Hamiltonian for the z-direction. Requiring that B is strong, we 
can approximate (j2~Tj) as 

„ m . h 2 k 2 z 1 (eB\ 2 2 . . 

H^\ sB = —^ + — — z 2 (106) 
2m 2m \ c J 

which is nothing but the one-dimensional harmonic oscillator of frequency u c , with the abbreviation 
(sB) means strong magnetic field. Therefore, the corresponding spectrum is given by 

r 



E n \sB = f^>c n + 



This modifies the eigenvalues as 

^n,a(kB)\ s B = ef' (k B + (jk s 

as well as the matrix element 

{n\k z \n) = ±i6 n>n f±i\ 



fc 2 k 2 c 
2m 



+ ftuj c I n + 



1 



' muj c 
~2h 



1 1 

n + - =F - 



(107) 



(108) 



(109) 



Note that, what makes difference with respect to the previous analysis is that the parameter Afc )Cr is 
absent. This of course will offer different simplifications and determine interesting results concerning 
the polarizations. 

To be coherent with our analysis, we consider the two first transverses modes and derive the whole 
spectrum. Doing this to get 



^,77=1,2(&b)|sB = £/ + (^b) - 6i, 2 (^b)|sB, £/,r;=3,4(feB)|sB = (k B ) + 62,1 IsB 

where different quantities changed to 

£i~(k B )\ B B = ^ {4 0) ( k B + k so) + e| 0) (k B - k so )} + ftw c 



h 2 k 2 so 

2m 



as well as 



and we have set 



£1 



h k so moj, 



sL! 



■;??, 



2h 



(110) 

(111) 

(112) 
(113) 



The corresponding Bloch spinors are given by 

^,fc B ,j7=l,4|sB = 



V 





1 



(114) 



/ 
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and the other sub-band read as 



:,fc s ,r/=2,3|sB 



/ \ 



V 



i 





(115) 



/ 



Is relevant to discuss the polarization symmetries in the present case to underline what makes 
differences with the previous analysis. On the other hand, the above analysis is sharing some commun 
features with that realized in |15| . Indeed, using the above tools to show that for r/ = 1,4 we have 



pWi 



N l,k B ,r,=l,A X 



£ l ( k B) ~ -7T T&i(&b)|sB 



(116) 



as well as for r\ = 2, 3 



P {Z) \ R 



Finally, we end up with 



N, 



-l 

kg,ri=2,3 



, X 1 



T^r- (k B ) + ^ ± £,i 2 (k B )\sB 
s IsB ^ 



p£L 1)4 (^)| s b 



-^ ) =2,3(-^)|sB. 



(117) 



(118) 



We reach the conclusion that the polarizations for strong magnetic field are conserved. This is an 
interesting results because it gives another way to talk about such symmetries, which is different from 
what has been developed in [15] . In fact, these polarizations is actually controlled by the external 
parameter B and can be adjusted by switching on B to diffrent values. 



9 Critical point analysis 

As we have noticed before in diagonalizing the Hamiltonian of the present system there is a critical 
point, i.e. Afc j(J = ±\J\, with a = ±1. It is a good task to study the effect of this point on the obtained 
results and underline what it will be new. To do this, we return to the Hamiltonian describing a 
system like harmonic oscillator where Afc j(7 is included and resort the corresponding spectrum to see 
what makes differences with respect to the standard case. 

Taking into account such critical point, we can establish a relation between the Bloch wave vector 
and the Rashba one. More precisely, we have 



mu> c mu> c 
k so \ a=1 = Ty-^ft- - k s, k so \ a =-x = ±J -— + k B - (119) 

It is convenient to write these relations as functions of the magnetic length l B , such as 

k so \<T=i = =F^=: k B , k so \ a= _i = ± — — + k B . (120) 

V 2l B V 2L B 

These relations are interesting in sense that one vector can be determined in terms of the other. 
Moreover, they are -B-dependent and clearly for a strong B one of them can be ignored where the 
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remaining vector will be controlled by B. In this limit for instance the area of the Hall droplet can be 
written as 

S = 2vr/|^ I (121) 
Now let us return to write the corresponding Hamiltonian, which is 

HP\ fI = huj c a J[ a. (122) 

This leads to end up with the spectrum 

/„f\n 

42,^)1 + /T = ^n, |n) = |0>, n = 0, 1,2 ■ ■ ■ . (123) 

The first remark one can underline is that actually the Hilbert space of such spectrum has an iso- 
lated point that correspond to a zero mode energy. This is not a suprising result because there are 
some systems who are behaving like this. In fact, electrons in graphene, which behave like massless 
Dirac fermions (Majorana fermions), in the presence of a magnetic field have such kind of zero mode 
energy [19] . More precisely, for one single particle in graphene can be described by the Dirac operator 



i 



S 



Ha = ^\- d o ) <124) 

where Vf ~ is the Fermi velocity, which will be set to one, and the many-body effects are neglected. 
The annihilation and creation operators, in the complex notation z = x + iy, are given by 

d B + d B . . 

d=2 a; + r- d, = - 2 8i + 2 3 - < 125) 

They verify the commutation relation 

[d,S]=2B. (126) 

Now let us look at its square, which is 

Clearly, the first matrix element is similar to that of (|122p . The corresponding wavefunctions should 
be written in an appropriate form. This is 

*m,n = ( ^ ) (128) 



•0, 



m— l.n 



where the eigenfunctions ijj m ,n are given by 



V> m ,„(z, z) = ( irV^ m ^! n L n fzz\ e _| Kgf m , n = 0, 1, 2 • • • . (129) 



^2"+ 1 7r(m + n)! V 2 

Their Landau levels take the form 

El(m)=Bm. (130) 
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Finally, we end up with the zero mode energy for such system, namely 



V>0,n 





(131) 



This is showing the similarity with the previous case that corresponds to the Hamiltonian at the 
critical point. 

Now let us return to derive the whole spectrum in the present case. Indeed, as an immediate 
consequence the eigenvalues of becomes 



tflML R = ef ] (k B + ak so ) 



h 2 k 2 

— h hio c n. 

Zm 



(132) 



The matrix element (1561) now reads as 



,,,,/, , . c / muj c ( 11 

{n\k z \n ) = ±td n ^ n i±i a —— I n + - T - 



(133) 



Focussing on the two first transverse modes, we show that different quantities forming the total 
spectrum take other forms. Indeed, we find 



hjjj r h 2 k 2 so 

2m 



as well as 



where £ at the critical point is given by 



- 1 - \ / o 



h 2 k. 



mui c 



m \ 2h 

These leading to end up with the energies for different sub-bands, namely 



(134) 



(135) 



(136) 



^=1,2 (&b) = tt{k B ) - ii{k B )\ /y, e ltV=3j 4(k B ) = ef{k B ) + £i(k B )\ n 



The corresponding Bloch spinors read as 

/ 



^,fc B ,r?=l,4 



si 



e;{k B ) + ^+ti(k B ) 



1 



(137) 



(138) 



as well as for other modes 



7 Z,fc s ,r)=2,3l ± R 



( 



s-(k B ) + ^-^(k B ) 

1 





(139) 



19 



The above materials can be used to analyze the polarization symmetries at the critical point. In 
fact, we show that for longitudinal and perpendicular cases the polarizations are zero. However it is 
finite along z-direction and this can be traduced by the obtained results. Indeed, for r\ = 1, • • • ,4 we 
obtain 

1 



V 2 V 2 



5 l± /i 



and for r/ = 2, 3 we have 

>(*) 



P 



N- 



l,r) l _i_ /i i,fcs ,77=2,314. /l 
V 2 V 2 



x < 1 



huj, 



1 2 



1 



n 2 



±^2 



(140) 



(141) 



These allow us to end up with the relation 



P 



S=i,4(Mi ± yr + -^;u, 3 (-^)i ± yT 



(2) 



(142) 



which means that there is no preservation of such symmetries as we have seen in the standard case. 
With this we end the analysis of four cases and we will see the last one. 



10 Periodic structure without magnetic field 

To complete the present analysis it is important to treat the case of the periodic structure without 
magnetic field. This will allow us to get new results and make contact with other proposals. In doing 
so, we simply return to the total Hamiltonian and derive the corresponding spectrum when B is zero. 

To start let us consider the Hamiltonian describing the system in the absence of B. This can be 
obtained from (fTT|) to end up with 

tap tat. 

H\ B=0 = - s JL{ ax k z -cj z k x ) + U{x). (143) 

2m m 

Clearly for arbitrary periodic potential U(x), the eigenvalues (|35j) and eigenstates (f36l) reduce to those 
of wave planes. Indeed, it is easy to notice 

These can be used to show that the matrix element takes the form 

(k z \k z \k' z ) = i5 kztK k z . (145) 

Combining all together to get the energies corresponding to different sub-bands. Note that, we have 
now only two possible values for 77, i.e. (1,2). Thus, we have 



£i, v =i,2(kB, k z ) = - ej 0) (k B + k so ) + e[ 0) (k B - k so ) 



h 2 kj h 2 k 



2m 2m 



±\i[sr(k B )) • mo) 



in 
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At this point, there are two limits those should be discussed to characterize the present situation. 
The first one is concerning the case where we assume that ks = and k z > 0. Returning to the 
spectrum to obtain 

,^ h 2 k 2 h 2 k 2 h 2 k k 

^=iA k B = 0, k M ) = e °> (k„) + - ^ ± (147) 
' ' ' 2m 2m m 

This tell us that the energy branch with rj = 2 has a minimum at k z = k so for all band I. These 

analytical results are in contrast to what was numerically predicted in |20j . 

The second limit can be fixed by requiring that U(x) = 0. In this case, one can show that the 

corresponding energies are 

h 2 / \ 
e„, v =i,2{k x , k z ) = — (k 2 x + k\ + 2ak so y / k 2 + k 2 ) . (148) 

In conclusion (|148j) is in agreement with the Rashba energy for the present system without periodic 
potential and magnetic field. 

11 Conclusion 

We have investigated the basic features of two-dimensional periodic structures with Rashba spin-orbit 
interaction in the presence of an external magnetic field B. This latter allowed us to end up with a 
confining potential along z-direction that has been used to deal with different issues. In particular, it 
has been served to split the corresponding Hamiltonian into two parts and one of them is obtained to 
behave like one-dimensional harmonic oscillator. This decomposition was usefull in sense that different 
spectrum's are obtained and leading to the total one. To make comparisons with other proposals, we 
have restricted ourselves to the two first transverse modes. These are used to derive the eigenvalues 
and the Bloch spinors for four sub-bands r] = 1, • • • ,4. Looking at the obtained spectrum we have 
shown that there is an hidden symmetry, which has been used to talk about the polarizations. 

After getting the eigenvalues and the corresponding eigenspinors, we have analyzed the polarization 
symmetries by evaluating their components in different directions. As consequence, the longitudinal 
and perpendicular are obtained to be nulls except that along z-direction where a finite values is found. 
Looking at the symmetry with respect to the Bloch vector, we have shown that the presence of the 
magnetic field breaks such symmetries in disagreement with the results obtained in [15J. As another 
task, we have investigated the dynamical spin by determining some relevant quantities those used to 
write down the appropriate spin Hall conductivity. 

Subsequentely, different limits have been studied. Indeed, we have analyzed the case when the 
Rashba coupling is null and this allowed us to end up with a Landau problem submitted to the 
periodicity constraint. Therefore, the spectrum is easily obtained by noting that the corresponding 
energies are degenerates and their degree of degeneracy can be determined in terms of the Bloch vector. 
Moreover, we have considered another case that is U(x) = 0, which lead to get the Landau problem 
amended with the Rashba coupling in the presence of the magnetic field. Solving the eigenvalue 
equation we have derived the spectrum that has been used to conclude that the polarizations are not 
preserved. 

Inspecting other limits, the strong magnetic field case has been considered. This allowed us to get 
different conclusions and derive interesting results. Indeed, after getting the corresponding spectrum, 
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we have tackled the symmetries problem. More precisely, we have shown that in the present situation 
the polarizations are conserved and therefore offered another way to talk about such symmetries. 
More importantly, such analysis allowed us to make contact with the proposal developed in [15] and in 
particular to show how one can recover the symmetry preservation's from the present work. Moreover, 
the obtained polarization's are magnetic field dependent and therefore it can be adjusted by switching 
such field to any value. 

On the other hand, we have studied in details one important issue that allowed us to make contact 
with other systems. In fact, the critical point Afc j(J = ±\J\ that appeared as a singularity in diagonal- 

izing the Hamiltonian has been analyzed. Resorting the corresponding spectrum we have ended 
up with a zero mode energy in the Hilbert space of H { z 1] . In fact, a comparison to the Dirac fermions 
in the presence of the magnetic field is established. 

Finally, to get more information about the present system we have investigated the case where B 
is switched off. This has been done by deriving the spectrum for r\ = 1,2, which are used to study 
two limits. Indeed, assuming that ks = where k z > 0, we have shown that the energy for r\ = 2 has 
a minimum at k z = k so for all band I, which is in agreement with that has been numerically predicted 
in [20]. Moreover, requiring that U(x) = 0, we have recovered the Rashba energy. 

Acknowledgment 

This work was completed during a visit of AJ to the Abdus Salam Centre for Theoretical Physics 
(Trieste, Italy) in the framework of junior associate scheme. He would like to acknowledge the financial 
support of the centre. 

References 

[1] For instance see R.E. Prange and S.M. Girvin (editors), "The Quantum Hall Effect", (Springer, 
New York 1990). 

[2] K. von Klitzing, G. Dorda and M. Pepper, Phys. Rev. Lett. 45 (1980) 1405. 
[3] D.C. Tsui, H.L. Stormer and A.C. Gossard, Phys. Rev. Lett. 48 (1982) 1559. 
[4] M.I. D'yakonov and V.I. Perel, JETP Lett. 13 (1971) 467. 

[5] Y. Kato, R.C. Myers, A.C. Gossard and D.D. Awschalom, Science 306 (2004) 1910. 

[6] J. Wunderlich, B. Kaestner, J. Sinova and T. Jungwirth, Phys. Rev. Lett. 94 (2005) 047204, 
cond-mat/0410295| 

[7] H-A. Engel, E.I. Rashba and B.I. Halperin, Theory of Spin Hall Effect in Semiconductors, 
cond- mat/0603306, 

[8] G. Feve, W.D. Oliver, M. Aranzana and Y. Yamamoto, Phys. Rev. B66 (2001) 
155328 fcond-mat/010802l| 

[9] M. Governale and U. Ziilicke, Phys. Rev. B66 (2002) 073311. 



22 



[10] V.Ya. Demikhovskii and A. A. Perov, Europhys. Lett. 76 (2006) 477, |cond-mat/0607216[ 

[11] B. Molnar, P. Vasilopoulos and F.M. Peeters, App. Phys. Lett. 85 (2004) 612; ibid Phys. Rev. 
B71 (2005) 125301; X.F. Wang, P. Vasilopoulos and F.M. Peeters, Phys. Rev. B72 (2005) 
075330. 

[12] V.Ya. Demikhovskii and A. A. Perov, Phys. Rev. B75 (2007) 205307. 

[13] D.V. Khomitsky, Phys. Rev. B76 (2008) 033404; ibid B77 (2008) 113313. 

[14] F. Correa, V. Jakubsky, L-M. Nieto and M.S. Plyushchay, Phys. Rev. Lett. 101 (2008) 030403, 
larXiv:0801. 16711 [hep-th]. 



[15] S. Smirnov, D. Bercioux and M. Grifoni, Europhys. Lett. 80 (2007) 27003, [arXiv:0805. 3830 [cond- 
mat.mes-hall]. 

[16] K. Arii, M. Koshino and T. ando, Phys. Rev. B76 (2007) 045311. 



[17] B.A. Bernervig and S-C. Zhang, Phys. Rev. Lett. 96 (2006) 106802, cond-mat/0504147 



[18] J. Schliemann, J.C. Egues and D. Loss, Phys. Rev. B67 (2003) 085302. 

[19] A. Jellal and M. Bellati, Fractional Quantum Hall States in Graphene, arXiv:0805.2388 [hep-th]. 

[20] V.Ya. Demikhovskii and D.V. Khomitsky, Pis'ma Zh. Eksp. Teor. Fiz., 83 (2006) 399. 

[21] J. Luo, H. Munekata, F.F. Fang and P.J. Stiles, Phys. Rev. B41 (1990) 7685. 



23 



